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Abstract 



We compute the heavy-heavy, light-light and light-heavy quarkonium spec- 
trum starting from a first principle Salpeter equation obtained in a preceding 
paper. We neglect spin-orbit structures and exclude from our treatment the 
light pseudoscalar states which in principle would require the use of the full 
Bethe-Salpeter equation due to the chiral symmetry breaking problem. For 
the rest we find an overall good agreement with the experimental data. In 
particular for the light-light case we find straight Regge trajectories with the 
right slope and intercepts. The strong coupling constant a s , the string ten- 
sion a occurring in the potential and the heavy quark masses are taken from 
the heavy quarkonium semirelativistic fit with only a small rearrangement. 
The light quark masses are set equal to baricentral value of the current quark 
masses as reported by the particle data group. For what concerns the light- 
light and the light-heavy systems the calculation is essentially parameter free. 
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In reference [T[ a Bethe-Salpeter like equation was derived which provides in principle a 
fully relativistic description of the quark-antiquark system. The derivation takes advantage 
of an appropriate path integral representation for a kind of 4-point gauge invariant Green 
function and proceeds entirely from first principles apart from the so called modified area 
low assumption. Such assumption consists simply in adding an area term to the perturbative 
value of the logarithm of the Wilson loop correlator, and, as well known, it is well founded 
both on lattice simulation and on theoretical considerations 



i\nW = i(]nW) peit + aSi 



(1) 



Unfortunately a direct treatment of the full B-S equation in its four dimensional form 
seems to be beyond the present possibilities. From it, however, by a standard method (which 
includes replacement of the full quark propagator by a free propagator and an instantaneous 
approximation of the kernel) a center of mass hamiltonian can be obtained in the form 



H = J m\ + k 2 + J ml + k 2 + V, 
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In eq. (@4J) the perturbative term has been evaluated only at the first order in the 
coupling constant a s , the indices 1 and 2 denote the quark and the antiquark, k and k' 
denote the final and the initial center of mass momentum of the quark, a is the string tension, 



k+k' . 



; Q = k' - k; g iC 



Wj+ttl'- 



= (5 hk — f h f k )q k is the transverse momentum, while 



a k are the usual Dirac matrices 7^7^, and a k the 4x4 Dirac spin matrices i/4e knm [j'j, 7™]. 

It is shown in [|| that by an 1/m expansion and an appropriate Foldy-Wouthuysen 
transformation one arrives from eqs. (|3|) and (|J) to the semirelativistic potential discussed 
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in 0. Such potential is made by a static part Vq, a velocity dependent part V v d and a spin 
dependent part V^d- The expression V s ^ is identical to that derived from the so called scalar 
confinement hypothesis while V wt ± is different. 

As shown in [0 and |§, when appropriate values are given to the quark masses rrij, to a s 
and to a, the semirelativistic potential Vq + Kd reproduces reasonably well the spin averaged 
multiplets in the heavy quarkonium case and it is phenomenologically favoured with respect 
to the scalar confinement expression. On the other hand it is shown in []5| that V s d gives the 
correct fine and hyperfine splittings if even the second order contributions in a s are included. 

In this paper we want to treat directly the hamiltonian (0-^). In this way we shall see 
that we not only reproduce the same results for heavy quarkonia as from the semirelativistic 
potential, but we also obtain the correct ground Regge trajectories (with the right slope and 
intercepts) for the systems involving light quarks alone and the known lowest states for the 
systems with a light quark and a heavy antiquark. (On the contrary we did not succeed in 
evaluating daughter trajectories due to computation difficulties.) As we shall explain more 
precisely later, we have achieved this goal by a very small rearrangement of the parameters 
used in |^ and adopting current masses as given by the particle data group || for the light 
quarks u, d and s. 

Notice that to treat light-light and light-heavy systems various models have been at- 
tempted in literature, based on some conjectural choice of a B-S kernel (see for instance [0) 
or directly of a local potential in eq. (0) (see e.g. || and references therein). However while 
such models give a reasonable spectrum, they usually fail in reproducing the observed Regge 
trajectories (in the sense that they do not result straight @, or, if they do, they have not 
the correct slope and intercepts) for the values of the string tension required by the heavy 
quarkonia fit. 

In this paper we shall neglect as a rule the spin dependent terms in eqs. (|3|) and (|]) 
and shall mainly pay attention to the velocity dependent part alone. Indeed the inclusion 
of the spin-orbit terms in the light quark case would make the entire treatment much more 
involved due to the impossibility of an expansion in 1/m and, as we shall see, it would not be 
particularly significant for our analysis. Later we shall return on S-wave hyperfine splitting 
which however would require the inclusion of the a? terms to give quantitative predictions. 
Finally in the light-light case we shall restrict our attention to triplet states alone. In fact for 
the light pseudoscalar mesons like tt, K, 77, 77', the approximation of the quark propagators 
by their free expressions as implied in the use of an hamiltonian is certainly inadequate, due 
to the complicate interplay existing in this case between quark propagator and B-S wave 
function which is related to the chiral symmetry breaking |J. 

Restricting our consideration to the velocity dependent part of the potential and splitting 
it in a perturbative and in a confinement part, V = V pcrt + V con f, we can write 
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In this way our hamiltonian becomes strictly related, but not identical, to that considered 



in ref . llC 



Let us begin to consider potential @ and (H) in various different limit situations. First 
let us take the static limit consisting in setting k = k' = in d§-[?p. In this way fl5|) and (^) 
become local and the hamiltonian can be written as 
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Less drastically let us assume <C q 2 (small angular momentum) and expand (0) in such 
quantity. Introducing the regularized Fourier transform 
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L 2 being the orbital angular momentum of the two particles (notice that classically q^ 
, ' 1 - > L ~ 1 one finds 
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On the contrary setting in eq. (0) q^ ~ q 2 (large angular momentum) one has 
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Finally if in eq. (0) we consider the extreme case mi = m 2 = 0, neglect VJp er t and simply 
set k = k' in the factor multiplying the singular term 
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we obtain the hamiltonian 
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In literature eq. @ has been used as such. In ref. ||, e.g. an overall fit of the low angular 
momentum meson spectrum has been obtained by setting a = 0.22 GeV 2 , a s = 0.323 and 
adding to the lefthand of eq. (||) an ad hoc flavour dependent constant C in the sector of the 
light quarks. As well known, however, for small mi and m2 eq. (H) would give asymptotic 



straight Regge trajectories with slope a' = 1/8<t, which for a - 
GeV -2 while the experimental value is about a' = 0.88 GeV -2 . 
On the contrary from eq. (|i~4|) we would obtain 



a: 
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0.22 GeV 2 gives a' = 0.57 



(15) 



which is identical to the Nambu-Goto string model. This equation gives the correct experi- 
mental value for o ~ 0.18 GeV 2 , which is the value used e.g. in ref. ]3|,|[] to fit the cc and the 
bb spectra (when the coupling with the decay channels is neglected). This result illustrates 
the advantage of considering the potential (§-0) with respect to a simply local potential of 
the type appearing in (||). 

Actually we do not succeed in directly diagonalizing the hamiltonian eqs. (^],|5|-[?J) , due 
to problems of numerical stability and computer time. Therefore we follow the following 
strategy. First we diagonalize the static hamiltonian H stat by the Rayleigh-Ritz variational 
method using the harmonic oscillator wavefunctions as basis ]n|,[nj. These wavefunctions 
in coordinate space take the form 
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where A is a scale factor and L n _l are the Laguerre polynomials [13[]. Then we write the 
eigenfunctions of H stat as 
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(19) 

Finally the eigenvalues of H are estimated as the minima of ( |19D in A. 

Due to eq. ( |T7| ) the evaluation of the kinetic term in eq. ( |T9"D is simply reduced to a one 
dimensional integral. 

For the perturbative term, by setting 
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(Pi and being respectively the first and second kind Legendre polynomials) a double 
integral. 

For what concerns (ip n i m \ V con f \ VWm}> if we had used the original expression (|,0), we 
should have ended in a five- dimensional integral which numerical evaluation is quite prob- 
lematic for a strongly oscillating function like ip n im{^, k). For this reason in place of @ and 
(0) we take in turn advantage of the approximate eqs. ( |TTD or (0). From eq. (|ll|) we obtain 
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( again a double integral) and from eq. ([12|) 
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or 



(a triple integral). The singularities occurring in the integrals in eqs. fl2"T|) and 
can be handled by the method explained in ref. Jl4| . 

We use eq. (|2~2"1) for the evaluation of the heavy quarkonia spectrum and for the light 
quarks S state, eq. (f23|) for the Regge trajectories. 

Notice that eq. (p2|) produces systematically larger masses. The difference being of the 
order of few MeV for the bb system, between 5 and 15 MeV for cc and progressively larger if 
light quarks are involved. The subsequent term in expansion (^) would contain a L 2 (L 2 + 2) 
factor and the convergence of the expansion soon becomes very slow as L increases. For this 
reason we belive eq. (^) to be preferred for light quarks already for the P states. 

Once that the spin averaged masses are determined we may evaluate the S-wave hyperfine 
splitting by the equation 
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this is obtained taking into consideration only the pure <T\ ■ cr^ term in eq. (Q). 

We have adopted the following parameters: a s = 0.363, o = 0.175 GeV 2 , m c = 1.40 
GeV, nib = 4.81 GeV, m s = 200 MeV, m u = 10 MeV; no ad hoc constant C has been added 
to the potential. The first four values have to be compared with those obtained from heavy 
quarkonium fits. E.g. when pair creation effects are neglected (and after renormalization of 
the masses by reabsorbing the constant C) the values used in ref. (table 2, first column) 
become a s = 0.363, a = 0.178 GeV 2 , m' c = m c + § = 1.397 GeV, m' b 



m b + § = 4.792 



GeV, (cf. also ref. fig. 4). On the contrary the light quark masses are taken as the 



baricentral values of the current masses: 

, current 



m 



current 



2 to 8 MeV, 



m 



current 



5 to 15 MeV, 

777,~""~" = 100 to 300 MeV as reported from the Particle Data Group ||. Notice also in this 
connection m c ^ TTent = 1.0 to 1.6 GeV, m£ urrcnt = 4.1 to 4.5 GeV. The small rearrangement in 
the values of a, m c and mj is required to obtain the exact slope for the Regge trajectories 
without modifying the ip and the T ground states. Notice that apart from that no attempt 
of optimizing the parameter is made. 

The results of the calculation for the cc, bb systems are reported in table | and compared 
with the spin averaged multiplets. An estimate for such values is taken from ref. f|. Also 
the unperturbed values obtained by H stat are reported, as it can be seen the shift for heavy 
quarks is of some tenth MeV (but such shift becomes more important for light quarks). 

As expected we have reasonable overall agreement with the data, the discrepancies being 
in part ascribable to pair creation effects ||. 

In Fig. [I], |3] the Regge trajectories are reported for the ground triplet states of the 
systems uu, ss, us corresponding to J = L+l, L, L — l. The agreement is again usually very 
good, particularly for the states J = L + 1, while some disagreement, particularly in the us 
case, can possibly be retraced in the neglecting of the spin-orbit terms. In table numerical 
results are reported for the uu system and compared with the data for the J = L + 1 state. 
The values in bracket for the IS and 2S states are obtained by adding 1/4 of the hyperfme 
splitting as given in table [jV[ Notice that the 2S state of table || is too high by about 100 
MeV, such discrepancy is however smaller than the width of p(1450) and u;(1420) and again 
could be ascribable to pair creation effects. Similar circumstances occur for the us and ss 
cases. 

The results for uc, ub, sc, and sb systems are reported in table |T| and compared with 
the experimental spin averaged masses using the theoretical splitting when the singlet state 
has not yet been observed. 

In table [TV] the hyperfme splittings as evaluated by eq. (|24J) are reported for the heavy- 
heavy and the heavy-light cases. In the light-light case a comparison with the splitting p-ir, 
ip-rj' and K*-K would have no meaning for the reason we have explained. Notice however 
that for the uu ground state eq. (|24| ) would give A = 221 MeV in good agreement with the 
uo-7] splitting (234 MeV). 

The fact that we do not have to use constituent masses in our calculation seem surprising 
at first sight. Notice however that, following |15 
hamiltonian eq. 
we find 



if we expand the kinetic part of the 
around the expectation value {ipniml^l^nim) (that we call simply (k 2 )) 
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M = m leff + m 2cfr + + — -3 — — + V + C (25) 

2mi eff 2m 2e ff 8mf efr 8m$ eS 



where at this order 



(k 2 ) 



m ieS = + (k 2 ) ^1 + - — j , i = i, 2 (26) 

and the quantity C can be defined by difference. 

Average values of m ue s and m se ff for the lowest states involving quarks u and s obtained 
evaluating (k 2 ) for the wave function eq. (|i~8|) are typically m uc s = 360 MeV, m se s = 410 
MeV, m ccS = 1.440 GeV, m bcS = 4.820 GeV, with C ~ 300 MeV. These are of the same 



order of the constituent masses used in semirelativistic computations [15]. 

In conclusion, starting from our first principle Salpeter equation, we have obtained an 
overall good reproduction of the spectrum of the mesons involving heavy and light quarks 
with the exception of the light pseudoscalar states. Since the parameter are practically 
completely specified by the heavy quarkonium spectrum and by high energy scattering, our 
calculation is essentially parameter free, for what concerns light-light and light-heavy quark 
systems. 

We would like to thank K. Maung Maung, L. Sorrillo and R. Chen for the useful discus- 
sions and N. Brambilla and A. Vairo for assistance on numerical problems. 
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TABLES 



TABLE I. Spin averaged energies for heavy quarkonium systems. Our theoretical results and 
experimental values. 



States 


static notpntial 
(GeV) 


tot a 1 "D OtPTl 1 1 a 1 


p"XT)priinpvitril valnps 


bb 








IS 


9.504 


9.446 


9.448 ± 0.005 


2S 


10.055 


10.014 


10.017 ± 0.005 


3S 


10.385 


10.347 


10.351 ± 0.005 


4S 


10.652 


10.615 


10.580 ± ? 


5S 


10.886 


10.849 


10.865 ± ? 


6S 


11.109 


11.068 


11.019 ± ? 


IP 


9.972 


9.953 


9.900 ± 0.001 


2P 


10.310 


10.288 


10.260 ± 0.001 


cc 
IS 
2S 
3S 
4S 
IP 
ID 
2D 


3.133 
3.693 
4.100 
4.441 
3.551 
3.833 
4.201 


3.065 
3.623 
4.030 
4.372 
3.508 
3.794 
4.158 


3.067 ± 0.002 
3.663 ± 0.005 
4.040 ± ? 
4.415 ± ? 
3.525 ± 0.001 
3.770 ± ? 
4.159 ± ? 
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TABLE II. Theoretical results and experimental data for uu system. For the S states we 
write the theoretical averaged mass and in brackets the mass for the triplet state obtained with 



the hyperfine splitting collected in table IV. 
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L.oZo 
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1 97c; _i_ pi nr\r 


i.oioi in u.uuui 


1 3 D 3 


1.721 


w 3 (1667) 0~(3~~) 


p 3 (1690) l+(3— ) 






1.667 ±0.004 


1.691 ±0.005 


1 3 F 4 


2.031 


/ 4 (2050) 0+(4++) 


a 4 (2040) l-(4++) 






2.044 ±0.011 


2.037 ±0.026 


1 3 G 5 


2.296 




p 5 (2350) l+(5— ) 








2.330 ± 0.035 


1 3 H 6 


2.531 


/ 6 (2510) 0+(6++) 


a 6 (2450) 1"(6 ++ ) 






2.510 ±0.030 


2.450 ±0.130 



TABLE III. Theoretical results for uc, ub, sc, sb systems (MeV). Experimental data are en- 
closed in brackets. 



State 


uc 


ub 


sc 


sb 


IS 


1973 (1973 ± 1) 


5326 (5313 ± 2) 


2080 (2076.4 ± 0.5) 


5418 (5404.6 ± 2.5) 


2S 


2600 (2623±?) a 


5906 (5897±?) a 






IP 


2442 (2438±?) b 


5777 (5825 ± 14) c 


2528 (2535.35 ± 0.34) 


5848 (5853 ± 15) 



a Obtained from preliminary Delphi data m{D*') = 2637 ± 8 MeV, m{B*') = 5906 ± 14 MeV 
TBI subtracting 1/4 theoretical hyperfine splitting reported in table |V[ 



^Estimated from m(D* 2 ) = 2459 ± 4 MeV, m(A) = 2427 ± 5 MeV. 
'From preliminary Delphi data |16| . 



TABLE IV. Theoretical results for qq hyperfine splitting (MeV). Experimental data are en- 
closed in brackets. 



State 


uc 


ub 


cc 


bb 


sc 


sb 


IS 


111 (141 ± 1) 


59 (46 ± 3) 


97 (117 ±2) 


102 


108 (144) 


60 (47 ± 4) 


2S 


59 


38 


59 (92 ± 5) 


42 


62 


40 
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FIGURE CAPTIONS 



Fig. 1 

Ground triplet uu Regge trajectories. Theoretical results (full line) compared with experi- 
mental data (circlet). Cross denote less established masses. 

Fig. 2 

Ground triplet ss Regge trajectories (with the same notations of fig. 1). 
Fig. 3 

Ground triplet us Regge trajectories (with the same notations of fig. 1). 
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